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Starting from the NRQCD Lagrangian the heavy quark-antiquark potential is written in terms of field strength 
insertions on a static Wilson loop. The relevant matching coefficients are given at the present status of knowledge. 
The short-range, perturbatively dominated, behaviour of the spin-dependent terms is discussed. 



1. NRQCD AND THE WILSON LOOP 
FORMALISM 

Heavy quark bound states provide an ex- 
tremely difficult but at the same time appealing 
system to test QCD. The difficulties are obvious. 
One is the mixing of different energy scales. This 
is a typical feature of any bound state problem 
in quantum field theory and makes tricky even a 
purely perturbative solution of it. An other con- 
ceptual difficulty is connected with the nonper- 
turbative nature of low-energy QCD. This sug- 
gests that nonpertubative contributions have to 
be taken into account in almost all QCD bound 
states. The reason why heavy mesons are appeal- 
ing is that the existence of an expansion parame- 
ter (the inverse of the mass m in the Lagrangian 
and the velocity v of the quark as a dynamical 
defined power counting parameter) makes possi- 
ble to handle the first difficulty and to keep under 
control the second one. The tool is provided by 
NRQCD §. This is an effective theory equivalent 
to QCD and obtained from QCD by integrating 
out the hard energy scale m. The Lagrangian 
comes from the original QCD Lagrangian via a 
Foldy-Wouthuysen transformation. The ultravi- 
olet regime of QCD (at energy scale m) is pertur- 
batively encoded order by order in the coupling 
constant as in the matching coefficients which ap- 
pear in front of the new operators of the effective 
theory. This ensures the equivalence between the 
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effective theory and the original one at a given or- 
der in 1/to and as- At order 1/m^ the NRQCD 
Lagrangian describing a bound state between a 
quark of mass mi and an antiquark of mass m2 
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This is the relevant Lagrangian in order to cal- 
culate the bound state observables up to order 
O («■*). A discussion of the operators appearing 
in (0) in terms of powers of the quark velocity 
can be found in 0,^. The coefficients C2"''', c"/\ 
... are evaluated at a matching scale fi for a par- 
ticle of mass mj. 

Nonperturbative contributions to the heavy 
meson observables can be evaluated directly from 
the Lagrangian (|l|) via lattice simulations [Q. 
Typically, since the hard degrees of freedom have 
been integrated out explicitly, the needed lattice 
cut-off n is expected to be larger (smaller in terms 
of energy) than the usual one with a clear reduc- 
tion in the computation time. Despite the ad- 



vantages, there are also some drawbacks in this 
method. In particular in this way we do not learn 
very much on our "analytic" knowledge on the 
QCD vacuum structure. Moreover computations 
on coarse lattices are not always under control. 
Therefore it is worthwhile to use the Lagrangian 
of NRQCD as a starting point and to work out 
the quark-antiquark interaction in the so-called 
Wilson loop formalism The advantage in do- 
ing so is that all the nonperturbative dynamics 
will be contained in gauge field averages of field 
strength insertions on a static Wilson loop. These 
can be very easily evaluated by means of some 
QCD vacuum model or via traditional lattice 
simulations 0| providing in this way a powerful 
method in order to discriminate between different 
models. The derivation of the quark-antiquark 
potential in the Wilson loop formalism from the 
NRQCD Lagrangian was first suggested in this 
context in Q and is discussed with details in 
Here we present only some results. The heavy 
quark-antiquark potential (assumed that it ex- 
ists) is given by 
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The "potentials" Vi, V2, ... are scale dependent 
gauge field averages of electric and magnetic field 



strength insertions on the static Wilson loop and 
are explicitly given in W is the gauge av- 

erage of the non-static Wilson loop. The expan- 
sion of it around the static Wilson loop gives the 
static potential Vq plus velocity (non-spin) de- 
pendent terms. S*^-') and L^-'^ are the spin and 
orbital angular momentum operators of the par- 
ticle j. The matching coefficients are defined as 
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= Cp^g ± Cp^g and d is the relevant contri- 



bution to the mixing coming from the four quark 
operators in Eq. (|l|) and will be given in the 
next section. Apart from the matching Eq. (H) is 
equivalent to the potential derived in In the 
next section we will give explicitly the matching 
coefficients and discuss briefly the relevance of the 
matching in order to have a short range consistent 
potential. 

2. MATCHING COEFFICIENTS 

Since for reparameterization invariance cs = 
2ci? — 1 H , all the spin dependent potentials given 
in Eq. (^ turn out to depend only on cf (if the 
mass of the particle is irrelevant we will omit to 
indicate it). This coefficient is known up to two 
loop in the anomalous dimension 
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where (3j are the usual /^-function coefficients, 
70 = 2Ca, 71 = 68C1/9 - 26CANf/9, ci = 
2{Ca + Cp), Nf is the number of flavors, Cp 
is the Casimir of the fundamental representation 
and Ca is the Casimir of the adjoint representa- 
tion. At the lattice scale used in the numerical 
values of this coefficient at the bottom and charm 
mass are cpimi,) ~ 1.06 x (1 -|- 0.15) = 1.22 and 
cpimc) ~ 1.27 X (1 4- 0.25) = 1.59 respectively. 
The corrections due to the one loop matching are 
relevant (15 % in the bottom case and 25 % in the 
charm case) and therefore of the same order of 
the next power in the velocity in the Lagrangian 
(|l|) (usually accepted values are (w^) ^ 0.07 and 
(v'c) - 0.24). ^ 

An evaluation of the coefficient cd associated 
with the Darwin term in the NRQCD Lagrangian 



is given in ||l0[| : 
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This corrects a previous wrong evaluation given in 
|8|. At the lattice scale used in |^ the numerical 
values of this coefficient at the bottom and charm 
mass are c_D(mb) ~ 0.76 and C£)(toc) — —0.08 re- 
spectively. As pointed out in , since the poten- 
tial AVq® manifests a 1/r behaviour, this term 
gives a flavor-dependent contribution to the cen- 
tral potential. However this contribution is sup- 
pressed in the bottom case by the bottom mass 
(see Eq. (||)) and in the charm case by the small- 
ness of the corresponding matching coefficient. 

Finally, the contributions coming from the 
four-fermion operators are usually suppressed ei- 
ther in as or in powers of the quark velocity 
V Nevertheless under RG transforma- 

tion the contribution to the spin-spin potential 
coming from the chromomagnetic operator in the 
NRQCD Lagrangian mixes with some of the lo- 
cal four quark operators. In order to take into 
account this mixing the delta contribution to the 
spin-spin potential has been added in Eq. (||) 
though it would be suppressed in as- The coeffi- 
cient d has been evaluated in 0]: 
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As noticed in ||8|] the presence of the match- 
ing coefficients in the expression for the potential 
makes possible the agreement in the short range 
region between the potential derived here with 
the traditional QCD one loop perturbative calcu- 
lation, e.g. in Let us focus on the spin-orbit 
terms Vi and V2. Comparing properly with [ pT| 
we get for V2 the perturbative contribution 
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where je is the Euler constant. This expression 
agrees very well with the lattice measurement of 



the same quantity shown in Fig. |l|. In the same 
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Figure 1. The spin dependent potential V2 as 
given by the lattice measurement of M . 



way we get for Vi the perturbative contribution 
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It is extremely interesting to compare the above 
expression with the short-range behaviour of the 
Vi potential as given by the lattice measurement 
shown in Fig. |^. Apart an overall shift propor- 
tional to the string tension and therefore of non- 
perturbative origin the agreement is very good. 
This is quite significant since the perturbative 
part of Vi is entirely due to loop corrections. As 
a consequence Vi is more sensitive than V2 to the 
matching scale /x. Notice that at very short dis- 
tances the function —V{, just because the log(/ir) 
term, is expected to become negative, but up to 
now no lattice data are available in this region. 
As a last comment we notice that due to the so- 
called Gromes relation -¥{ = VJ,' |l| a V5 
potential (in the notation of ) emerges also in 
Eq. (^ by collecting the contributions coming 
from the fourth and fifth line. The perturbative 
expression we get agrees with that one given in 
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Figure 2. The spin dependent potential —V{ as 
given by the lattice measurement of 0. 



3. CONCLUSIONS 

In the framework of NRQCD and at the present 
status of the matching we have given the expres- 
sion for the heavy quark potential in terms of field 
strength insertions on a static Wilson loop. This 
has the advantage that traditional lattice calcu- 
lations can be used in order to evaluate nonper- 
turbative contributions. Moreover in this way a 
comparison between different QCD vacuum mod- 
els can be performed directly in terms of Wilson 
loop expectation values. This approach has been 
developed with some extent in [|. Here we have 
emphasized the role played by the matching co- 
efficients in order to make consistent the short 
range behaviour of the potential that we obtain 
with the usual scattering matrix derived poten- 
tial. We noticed that present lattice data are sen- 
sitive to one loop corrections and to the matching 
scale. 

As a conclusion, let us mention two open prob- 
lems. In order to have a 10 % accuracy in the 
quarkonium spin splitting it is necessary to add 
to the Lagrangian (|l|) higher order operators ||]. 
The inclusion (if possible) of such operators in an 
expression like Eq. (||) is still to do. Moreover 
in order to obtain Eq. we have implicitly as- 
sumed the existence of a potential. Non-potential 
terms surely exist in perturbative QCD. How to 
treat it in a system affected by nonperturbative 



physics is still unclear. Interesting developments 
could come from a promising approach recently 
proposed in IT^ . 
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